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Abstract 

Available experimental data in respect to a possible formation of a strongly interacting 
mixed phase are outlined. A physical program is formulated for new facilities being opend in 
Dubna for acceleration of heavy ions with an energy up to 5 AGeV. 



Over the last 25 years a lot of efforts have been made to search for new states of 

\o : 

strongly interacting matter under extreme conditions of high temperature and/or baryon 



density, as predicted by Quantum Chromodynamics (QCD). These states are relevant 
to understanding the evolution of the early Universe after Big Bang, the formation of 
neutron stars, and the physics of heavy-ion collisions. The latter is of great importance 



since it opens a way to reproduce these extreme conditions in the Earth laboratory. This 
explains a permanent trend of leading world research centers to construct new heavy ion 
accelerators for even higher colliding energy. 

Looking at the list of heavy-ion accelerators one can see that after the first experiments 
at the Dubna Synchrophasotron, heavy-ion physics successfully developed at Bevalac 
(Berkley) with the bombarding energy to Ei a b ~ 2 AGeV, AGS (Brookhaven) Ej a & ~ 11 
AGeV, and SPS (CERN) Ei a i ~ 160 AGeV. The first two machines are closed now. The 
nuclear physical programs at SPS as well as at SIS (GSI, Darmstadt, ~ 1 AGeV) 
are practically completed. The new relativistic heavy-ion collider (RHIC, Brookhaven) is 
intensively working in the ultrarelativistic energy range ^s^n ~ 200 GeV to searching for 
signals of the quark-gluon plasma formation. In this respect, many hopes are related to 
the Large Hadron Collider (LHC, CERN) which will start to operate in the TeV region in 



2 



two-three years. The low-energy scanning program at SPS (NA49 Collaboration) revealed 
an interesting structure in the energy dependence of some observables at Ei ab ~ 20 — 30 
AGeV, which can be associated with the exit of an excited system into a decofmement 
state. This fact essentially stimulated a new large project FAIR GSI (Darmstadt) for 
studying compressed baryonic matter in a large energy range of Ei a \> = 10 — 30 AGeV 
which should come into operation after 2013 year |l|. The general scheme how available 
accelerators are spread over the world is presented in Fig^ 



Heavy Ion Accelerators 




FIG. 1: The present and future accelerators of relativistic nuclear beams. For each facility the 
name of the laboratory, the name of the accelerator and the energy range (center-of-mass energy 
per nucleon-nucleon pair ^snn in GeV) are given 

On the other hand, in JINR there is a modern superconducting accelerator, Nuclotron, 
which has not realized its planned parameters yet. The Veksler and Baldin Laboratory of 
High Energy has certain experimental facilities and large experience to work with heavy 
ions. This study may actively be supported by theoretical investigations of the Bogoliubov 
Laboratory of Theoretical Physics. The paper aims to give an overview of experimental 
results on manifestation of some peculiarities of strongly interacting QCD matter, ob- 
served in relativistic heavy ion collisions, and to present arguments that acceleration of 
heavy ions like Au at the Nuclotron up to the maximal planned energy Eiab = 5 AGeV, 
allows one to study properties of hot and dense QCD matter to be competitive at the 
world level. 

A convenient way to present a variety of possible states of strongly interacting matter is 
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FIG. 2: Artist's view of the phase diagram 
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a phase diagram in terms of temperature T and baryon chemical potential /i# (or baryon 
density ub), as shown in Fig|21 This schematic picture shows in which region of the 
diagram the given phase is realized and which colliding energies are needed to populate 
this region. In particular, the boundary of the deconfinement and/or chiral symmetry 
restoration transition may be reached even below bombarding energies planned in the 
FAIR GSI project but the nuclear matter compression should be high enough. 
This point is illustrated in more detail in FigH3 
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FIG. 3: Dynamical trajectories for central (b = 2fm) Au + Au collisions in T — hb (left panel) 
and T — fiB (right panel) plane for various bombarding energies calculated within the relativistic 
3-fluid hydrodynamics 3]. Numbers near the trajectories are the evolution time moment. Phase 



boundaries are estimated in a two-phase bag model. 



As is seen, a system, formed in a high energy collision, is fast heated and compressed 



4 



and then starts to expand slowly reaching the freeze-out point which defines observable 
hadron quantities. Indeed, at the maximal achievable Nuclotron energy = 5 AGeV 
the system "looks" into the mixed phase for a short time (the left part of FigOJ). However, 
uncertainties of these calculations are still large. If heavy masses for u, d quarks are 
assumed, the phase boundary is shifted towards higher (the right part in Fig|3J). On 
the other hand, these dynamical trajectories were calculated for a pure hadronic gas 
equation of state, and the presence of a phase transition may noticeably change them. In 
addition, near the phase transition the strongly interacting QCD system behaves like a 
liquid rather than a gas, as was clarified recently at small \x& from both quark 4| and 
hadronic 5] side. As to high fi B values, it is a completely open question. One should 
note also that as follows from lattice QCD calculations for /ig « the deconfinement 
temperature practically coincides with the transition temperature for chiral restoration 
symmetry while for baryon rich matter it is still an open question. 
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FIG. 4: Dynamical trajectory in the e — n#-plane for central Au+Au collisions calculated with 
two equations of state: pure hadronic (solid line) and with first order phase transition (dashed). 
Spatial averaging is done over the cube with the 4 fm sides and Lorentz contracted in the 



longitudinal direction 



. The shaded regions correspond to the mixed phase (upper one) and 



the non-reachable domain with the boundary condition T = 0, respectively. 



In FigEJ dynamical trajectories in the e — rig plane for the top Nuclotron energy are 
given for two equations of state, without and with the first order phase transition j^, 
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The shaded band corresponds to the mixed phase of more restrictive EoS with heavy 
quark masses, as shown in the right panel of FigEl Nevertheless, both trajectories, being 
close to each other, spend some time in the mixed phase. The main difference between 
the results presented in Fig|3] and Fig0] comes from the different averaging procedures 
used: in the first case, thermodynamic quantities were averaged over the whole volume 
of the interacting system, while in the second case, it was carried out only over a cube of 
4 fm sides placed at the origin and Lorentz contracted along the colliding axis. Thus, for 
central collisions at the 5 AGeV Nuclotron energy even if an average state of the whole 
strongly interacting system does not approach the mixed phase, an essential part of the 
system volume will have a certain time in this mixed phase. An experimental consequence 
is that an expected observable signal for reaching the mixed phase should be rather weak. 

It is hard to believe that some irregular structure like that at Ei a b ~ 30 AGeV []| can 
manifest itself at the Nuclotron energy. So the global observables are expected to be quite 
smooth with energy. Indeed, as is seen from the results presented in FigsEl and El covering 
the region of ~ 2-10 AGeV. The shape of rapidity spectra for newly produced mesons 
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FIG. 5: The tt + ,tt~ and proton rapidity distributions in central Au+Au collisions at bombard- 
ing energies 2, 4, 6, and 8 AGeV 8]. Solid lines are drawn along experimental points. For 
comparison, Gaussian distributions are given by dashed curves. 



is bell-like but cannot be described by a single Gaussian due to flow effects. Note that 
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a number of 7r + is not equal to that of tt~ though their difference essentially decreases 
with the bombarding energy. So the isotopic degree of freedom should properly be taken 
into account under theoretical consideration in the Nuclotron energy range. As to proton 
spectra, their shape is yet wider, ranging from the target to projectile rapidity. At E [ab ~ 
10 AGeV it is getting flat at the middle rapidity and may be described by a superposition 
of two Gaussian to SPS energies |3( and of three Gaussian at the RHIC energies. The 
shape of baryon rapidity distributions strongly varies with the impact parameter (see 
FiglHJ) taking U-shape for peripheral collisions. The noted trends are continued with 
increase in energy and governed mainly by the stopping power of colliding matter. To 
see any effect of phase transformation in a global rapidity distribution, a phase like that 
should be created in a large volume and should live for a long time. Possibly, it is seen 
as appearance of the third midrapidity source at the RHIC energies [in| . 
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FIG. 6: Proton rapidity distributions for different centrality of Au+Au collisions at 6, 8, and 
10.8 AGeV Q|. Filled symbols are measured points, empty ones are their reflection with respect 
to the y = y cm line. Numbers under distributions give the centrality percentage. 



Thus, sensitivity of global observables to possible phase transitions is expected to be 
weak, but it might not be the case for more delicate characteristics. In any case, due 
to the proximity of the phase diagram region under discussion to the confinement tran- 
sition and chiral symmetry restoration, some precursory phenomena cannot be excluded 
at a bombarding energy of about 5 AGeV, which opens a new perspective for physical 
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investigations at the Dubna Nuclotron. 

Below some arguments in favor of this statement are given: 

• Pro per ties of hadrons are expected to change in hot and/or dense baryon mat- 



ter 



2 



12j. This change concerns hadronic masses and widths, first of all for the 
cr-meson as the chiral partner of pions which characterizes a degree of chiral sym- 
metry violation and can serve as a "signal" of its restoration as well as the mixed 
phase formation. Rare decays in matter of vector mesons (particularly p and uj) are 
also very attractive. 
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FIG. 7: Left panel: e + e invariant spectra from central Pb+Au (40 GeV) a collisions [Lj|. Thin 
solid and dotted lines are hadronic cocktail and the calculated results for free p mesons, respec- 



tively. Appropriate thick solid and dash-dotted lines are calculated in the Rapp-Wambach 
and Brown-Rho scenarios. Contributions of different channels are shown as well. 



Q 



Right panel: Invariant mass distribution of dimuons from semi-central In+In collisions at the 



beam energy 158 j4GeV. Experimental points are from [13]. Solid and dashed curves are calcu- 



□ 



lated [16] in the dropping mass scenario using the p-mass modification factors as density and 
temperature-density dependent, respectively. Dash-dotted curve neglects any in- medium modi- 
fication. Dotted line indicates the hydrodynamically calculated p-meson decay at the freeze-out. 



The presence of in-medium modification of p-mesons has been proved in the CERES 
experiments, see an example in Fig|7| The observed essential enhancement of low- 
mass (0.2 < M < 0.7) lepton pairs, as compared to free hadron decays, is due to 
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the influence of hot and dense nuclear matter on properties of the p-meson spec- 
tral function. Both the Brown-Rho scaling hypothesis assuming a dropping p 
mass and a strong broadening, as found in the many-body approach by Rapp and 
Wambach result in reasonable agreement with experiment. However, a poor 
resolution in the dielectron mass prevents the discrimination of different physical 
scenarios of this effect in the CERES experiments. The situation is noticeably bet- 
ter in the recent NA60 experiment with muon pairs lfij . As is seen from the right 
panel of FigEl the quality of dimuon data is much higher allowing a good resolution 
of the spectral function. The depicted comparison with the dropping mass scenario 
shows a strong influence of the final result on the hypothesis used for in-medium 
dependence of the mass modification factor Qj. However, it seems to be impossi- 
ble to describe simultaneously dielectron and dimuon experimental data within the 
dropping mass scenario under the same assumptions. Similar data are absent for 
heavy ions in the Nuclotron energy range. 

As to in-medium a-meson decay, some indications were obtained in reactions in- 
duced by pions and photons III bja, Ilia] . In FigJHl the relative pion pair abundance 



Cff , defined as 



yA a A (M nw ) , a N (M n 



U rp (J rp 

is presented for different isotopic states of observed pions versus the invariant pion 



pair mass M nn . A sizable enhancement is observed at low M nn in the case of 

n 

7r + — ► 7r 7r reaction but not for the n + — ► tt + tt + case |17I |. This enhancement may 
be related to a possible tttt scattering in a nucleus via formation of a virtual scalar 
o meson which is forbidden for the ir + 7c + state. The effect is getting stronger for 
heavier nuclei. Recently, using the ELSA tagged facilities in Bonn the in-medium 
modification of the cu-meson has first been observed in the reaction r y+A — > uj+X —* 
, + lt lQ.A smal , shift in the u mass wa S detected. Note, however, that ho 
it- and 7-induced reactions we deal with comparatively low baryon density states, 
Ub ~ (1 — 2)n . There are no similar experiments with heavy ions. 

Nevertheless, there are theoretical proposals to probe chiral symmetry restoration 



in the vicinity of the phase transition boundary. In particular, it was shown 20L 12 1 1 



that a two-photon decay of the cx-meson formed as an intermediate state in ttit 
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FIG. 8: Invariant mass distribution of pion pairs from 7r-nucleus interactions 
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scattering may be a very attractive signal. As depicted in FigEJ at temperature in 
the vicinity of the phase transition, when m a ~ 2m n , there is an anomalous peak in 
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FIG. 9: Invariant mass spectra of 2j at \ib = 80 MeV and different temperatures 

invariant mass spectra of 7 pairs which may serve as a signal of the phase transition 
and formation of a mixed phase. Certainly, there is a huge combinatorial back- 
ground due to 7r° — > 77 decays, but the Nuclotron energy is expected to have some 
advantage against higher energy accelerators because the contribution of deconfined 
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quarks-gluons will be negligible. 

This effect may be observed in the e + e~ decay channel as well (3]. 

Electromagnetic probes discussed above carry out information concerning the whole 
evolution of colliding nuclei and states which are realized. The bulk of produced 
hadrons is related to the freeze-out point where information on the interaction 
dynamics has essentially been erased. So the expected behavior of global hadronic 
observables is smooth. However, some peculiarities of delicate hadron characteristics 
may be found and their hints are available even now. 
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FIG. 10: Excitation function of the proton elliptic flow [2j |. The shaded band corresponds to 
the Nucltron energy. 



In FigHni the beam energy dependence of the elliptic flow coefficient v 2 is presented 
for protons in the midrapidity range from noncentral heavy- ion collisions j^J. The 
elliptic flow characterizes the angular anisotropy in the transverse momentum plane 

as 

dN/dcf) ~ [1 + 2v! cos(0) + 2v 2 cos(20)] , 

where the <fi angle is measured from the reaction plane. As is seen, just near the 
Nuclotron energy Ei a b ~ 5 AGeV the coefficient t>2 changes its sign and the trans- 
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verse flow evolves from the out-of-plane towards in-plane emission. This fact may be 
treated as softening of the equation of state to be a precursor of the phase transition. 
The elliptic flow shows an essentially linear dependence on the impact parameter 
with a negative slope at Ei a h = 2 AGeV, with a positive slope at 6 AGeV and with a 
near zero slope at 4 AGeV [24]. This dependence serves as an important constraint 
to high-density behavior of nuclear matter for discrimination of various equations 
of state. 



Strangeness enhancement is an intriguing point of physics of heavy ion collisions, 
being one of the first proposed signals of quark-gluon plasma formation. An im- 
portant experimental finding is the observation of some structure ("horn") in the 
energy dependence of reduced strangeness multiplicity at Ei a b ~ 30 AGeV, pre- 
dicted in |25j as a signal that the formed excited system came into a deconfinement 
phase. As an example, in Fig^^the K + /ir + average multiplicity ratio is displayed 
as a function of the colliding energy. The "horn" structure is well visible. Such a 
structure is absent in free nucleon-nucleon collisions. 
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FIG. 11: Colliding energy dependence of the reduced strangeness abundance of K + mesons (2 
The shaded band corresponds to the Nuclotron energy range. 



The same experimental data (besides two recent points measured at E^ =20 and 
30 AGeV) are plotted in FigEH as a function of the bombarding energy. It is 
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FIG. 12: Beam energy dependence of the K + /tt + ratio. Curves are the results of different 
transport calculations (the figure is taken from 27]. The shaded band corresponds to the 
Nuclotron energy range 

of great interest that these ratios for strange hadrons are not explained by the 
modern transport theory (UrQMD, HSD models). While the mean pion and kaon 
multiplicities are well reproduced at the SIS and SPS energies, the above-mentioned 
models essentially underestimate the K/tt ratio in the AGS energy domain. It is 
remarkable that the divergence between the transport theory and experiment starts 
just at the Nuclotron energy. This increases interest in future experiments at the 
Nuclotron. 

• A strange particle ratio is not the only attractive observable. Let us look at the 
Lorentz invariant quantity, the transverse momentum distribution, presented in 
FigCni 

It is a global characteristic and no peculiarities are expected. Indeed, experimental 
spectra are reasonably well described within a simple blast-wave model [3] where 
spectra at freeze-out are parameterized as follows: 



dN 



T ^ ,rriT cosh p. sinhp, 
m T K l ( ) /„( 



dy rriTdrriT ' T ' T 

with the boost angle p = tanh -1 (3t and values of the freeze-out temperature T and 
transverse velocity (3t are given in FigH21 However, the behavior of the average 
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FIG. 13: The transverse mass spectra at midra pidi ty for different hadrons produced in central 
Pb+Pb collisions at 20 AGeV and 30 AGeV 26 1. Solid lines are results of the blast-wave 



parameterization 28] with the parameters given in the figure. 
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FIG. 14: The energy dependence of < tut > —nM) f° r pions, kaons, and protons at midrapidity 



for the most central Pb+Pb/Au+Au collisions 
energy. 



. Shaded bands correspond to the Nuclotron 



transverse mass, < ttit > —tuq, versus colliding energy y/s NN , is not trivial. As 
follows from Fig|T4l a remarkable change in the energy dependence around a beam 
energy of ~30 AGeV is clearly visible for pions and kaons exhibiting some kind of 
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FIG. 15: Variation of < my > with produced charged particles per unit rapidity at midrapidity 
for central collisions corresponding to the energy range from AGS to RHIC . The error bars 
reflect both the systematic and statistical errors in obtaining T e ff. 

plateau. While (< > —mo) rises steeply in the AGS energy range, this rise is 
much weaker from low SPS energies until RHIC energies where it starts again to rise. 
To a lesser extent this change is also seen for protons. One should emphasize that the 
beginning of the plateau is well correlated with the "horn" position. Measurements 
at the Nucloron may specify a pre-plateau behavior, in particular, for kaons and 
protons. 

Among many signals of the formation of Quark Gluon Plasma (QGP), one of the 
earliest is based on the relation of the thermodynamical variables, temperature, and 
entropy to the average transverse momentum and multiplicity, respectively, as was 
originally proposed by Van Hove 2t|. It was argued that a plateau in the transverse 
momentum beyond a certain value of multiplicity would indicate the onset of the 
formation of a mixed phase of QGP and hadrons, similar to the plateau observed in 
the variation of temperature with entropy in a first order phase transition scenario. 
In FiglTol the variation of < % > with charged multiplicity produced per unit 
rapidity dN/ dy around the midrapidty [3] is depicted for pions, kaons, and protons 
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FIG. 16: The energy dependence of the inverse slope parameter of kaons [3]. The shaded band 
corresponds to the Nuclotron energy. 

at AGS, SPS and RHIC energies spanning the range from 2 AGeV to 200 AGeV. 
In this analysis the experimental data on transverse mass spectra were parameter- 



ized as 



dN 



Cexp (-m T /T eff ) 



1) 



rriT drriT 

with the inverse slope parameter treated as an effective temperature T e ff. The 
average transverse mass of hadrons is related with T e f f : 



< m T > = T e ff + m + 



■eff) 



m + T eff 

From the results shown in Fig^Jone observes an increase in < > with a particle 
density dN c h/dy for AGS energies followed by a plateau for charge multiplicities 
corresponding to SPS energies for all the particle types considered: pions, kaons and 
protons. This may hint at the mixed phase, possible co-existence of the quark and 
hadron phases. For charge multiplicities corresponding to higher RHIC energies, 
the < rriT > shows an increasing trend indicating the possibility of a pure QGP 
formation. 

Using eq.Q, one may directly analyze data in terms of effective temperature. As 
is shown in FigEH the inverse slope parameter of kaons increases in the AGS and 
RHIC energy domains but it stays constant at SPS energies in natural agreement 
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FIG. 17: The colliding energy dependence of the mean multiplicities of various hadrons emitted 
in central Pb+Pb/Au+Au collisions 
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with the particle ratio results presented above. This feature, which is not observed 
in p+p interactions, might be attributed to the latent heat of a phase transition 25] 
and it is in fact consistent with hydrodynamic model calculations assuming a first 
order phase transition 31 1. 

As mentioned above, global multiplicities are expected to be quite smooth over 
a large range of energy. FigEI shows the mean multiplicities of hadrons emitted 
in central Pb+Pb/Au+Au collisions as a function of energy. The results are well 
described within the statistical hadron gas model 32 |. Note that a number of 
charged pions (n n + + n w - ) ~ 200 and kaons K + ~ per central event at the top 
Nuclotron energy. 

However one should add that peculiarities at the "horn" energy in strange particle 
ratios and the effective temperature dependence are also observed in average pion 
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FIG. 18: Pion multiplicity per participating nucleon for nucleus-nucleus (symbol) and nucleon 



33]. 



nucleon collisions (solid line) as a function of available energy in nucleon-nucleon collisions 
The shaded band corresponds to the Nuclotron energy. 

number per nucleon-nucleon interaction Such an example is given in Figll8l 

The reduced pion multiplicity < n n > /A part for nucleus-nucleus collisions is below 
that for elementary NN interactions at moderate energies and exceeds appropriate 
NN values in a relativistic domain. According to 0,0]) the equality of the nuclear 
and elementary < n n > /A part ratios at the "horn" energy is treated as an argument 
in favor of the onset of deconfinement. However, this proximity of nuclear and 
elementary reduced pion multiplicities starts at the Dubna Nuclotron energy and 
could naturally be explained by the role of the A isobar in the pion absorbtion. To 
clarify this or alternative interaction mechanisms, new experimental data for heavy 
ion collisions with scanning over the Nuclotron energy range are needed. 

• Existing experimental data indicate a large degree of equilibration of nuclear matter 
from high-energy heavy- ion collisions at freeze-out (see review- article 0])- The 
results presented in FigUm demonstrate a high quality description of relative hadron 
yields within the equilibrium statistical model obtained by means of fitting the 
temperature and baryon chemical potential. 

The mechanisms of equilibration are still unclear. In particular, a large difference 
between dynamical and statistical models is observed for heavy baryons (e.g. for A, 
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FIG. 19: Yield ratios with the best fit at midrapidity for the top AGS beam energy 10.7 AGeV 



351 ] . The diamonds 



(left panel) and for the RHIC colliding energy ^/snn =130 GeV (right panel) 
represent ratios of yields integrated over 4-7T. Note the scaling factor 100 for the ratios p/p and 
A/A in the left figure. The extracted freeze-out temperature and baryon chemical potential are 
displayed in the figures. 

multi-strange hyperons, H and Q~ and resonances) measured at high energies. The 
role of multiparticle interactions is also obscure. The corresponding data for the 
Nuclotron energies are extremely poor. Future Nuclotron measurements could give 
a strong input for a further development of dynamical and statistical approaches 
and allow one to disentangle these two types of models. 

• To a certain extent, the study of fluctuations in relativistic strongly interacting 
matter may help with solving the equilibration problem mentioned above. Experi- 
mental data on event-by-event fluctuations (e.g., fluctuations of particle multiplicity, 
electric, baryon and strangeness charges) in nuclear collisions give a unique possi- 
bility to test recent ideas on the origin of fluctuations in relativistic interacting 
systems Among them suppression of event-by-event fluctuations of elec- 

tric charge was predicted as a consequence of deconfinement. Estimates of the 
magnitude of the charge fluctuations indicate that they are much smaller in a quark- 
gluon plasma than in a hadron gas. Thus, naively, a decrease of the fluctuations 
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is expected when the collision energy crosses the threshold for the deconfinement 
phase transition. However, this prediction is derived under assumptions that ini- 
tial fluctuations survive through hadronization and that their relaxation times in 
ladronic matter are significantly longer than in the hadronic stage of the collision 
36. 38]. The magnitude of the measured charge fluctuations depends not only on 
the unit of electric charge q carried by degrees of freedom of the system (integer 
charge of hadrons or fractional one of quarks), but also on trivial effects, which may 
obscure the physics of interest (the fluctuations in the event multiplicity, caused by 
the variation of the impact parameter, and changes in the mean multiplicity due to 
changes of the acceptance). In our example, A$ q is used as a measure of charge fluc- 
tuation which is constructed from the well established measure <3> g of event-by-event 
fluctuations, defined as: 



$ 9 = v / < Z 2 > I < N > - Vz 2 (2) 

where: z = q — q, Z = J2iLi(li~ m these equations N is the number of particles 
of the given event within the acceptance, and over-line and < ... > denote averaging 
over a single particle inclusive distribution and over all events, respectively. 

For a scenario, in which particles are correlated only by global charge conservation 



(GCC), the value of § q is given by $ qi GCC = yl-P - 1 where P =< N ch > 
I < N ch > tot with < N ch > and < N ch > tot being the mean charged multiplicity 
in the detector acceptance and in full phase space (excluding spectator nucleons), 
respectively. In order to remove the sensitivity to GCC, the measure A$ 9 is defined 
as the difference 

A$ 9 = % - $ g ,GCC • 

By construction, the value of A$ g is zero if the particles are correlated by global 
charge conservation only. It is negative in the case of an additional correlation 
between positively and negatively charged particles, and it is positive if the positive 
and negative particles are anticorrelated 38]. 

As is seen from the results presented in FigJ2Hl the measured values of AQ q vary 
between and -0.05. They are significantly larger than those expected for QGP 
fluctuations —0.5 < A$ g < —0.15. A weak decrease of A$ g with increasing energy 
is suggested by the presented data. 
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FIG. 20: The colliding energy dependence of A$„ for a small (5y = 1.2 and large <5y = 3 rapidity 



interval for central Pb + P6 collisions 



3j|. Shaded bands correspond to the Nuclotron energies. 



Measurements require tracking detectors of large acceptance and precise control of 
collision centrality on event-by-event basis. Up to now only results on very limited 
acceptance at high energies are available, thus new measurements at the Nuclotron 
energy are of particular importance. From the experimental point of view the Nu- 
clotron energy range seems to be ideal for these measurements. This is because 
moderate particle multiplicity and their relatively broad angular distribution sim- 
plify an efficient detection of all produced charged particles. 



An investigation of the narrow Hanbury-Brown-Twiss correlations [40j joins this 
problem. Momentum correlations between mesons have been widely used to study 
the space-time extent of the emitting source in nucleus-nucleus collisions. For iden- 
tical bosons the symmetry requirement of the BoseEinstein statistics results in an 
enhanced production of bosons with a small momentum difference. Kopylov and 
Podgoretsky j^lj first noticed the deep analogy of this effect with Hanbury-Brown 
and Twiss spacetime correlations of the classical electromagnetic fields used in as- 
tronomy for interferometric measurement of star angular radii and developed the 
basic methods of momentum interferometry in particle and nuclear collisions. This 
technique represents an important tool for investigating the underlying reaction 
mechanism as far as it influences the decoupling configuration, including an unusu- 
ally long decoupling time suggested as one of the signatures of quark-gluon plasma 



formation 



421. 
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The correlation function C2(pi,P2) is defined as the ratio of the probability of ob- 
serving particles with four-momenta p\ and P2 simultaneously in one event divided 
by the probability for pairs of independent particles with the same single parti- 
cle phase space distribution. In Gaussian approximation with the Bertsch-Pratt 



parameterization 



42j the correlation function may be reduced to 



C^out, Qside, Qlong) — 1 + A exp [— (j^R^t ~ QsideRlide ~ ifongR] 



ong long] ' 



where qi ong is the longitudinal momentum component, and (q ou t, Qside) are the com- 
ponents of the transverse 3-momentum difference in an outward direction and per- 
pendicular to it. The parameter A characterizes the degree of incoherence. 
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FIG. 21: Bertsch-Pratt parameters as a function of colliding energy for 7r _ 7r~ pairs from central 
Pb+Pb/Au+Au collisions [^J. The presented data are near midrapidity and for comparable 
transverse momentum bins from nine different experiments. Systematic errors are not shown. 
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The parameters of the HBT analysis are presented in FigJ^ as a function of the 
colliding energy. Although there is some disagreement among experiments, the 
data do not appear to exhibit any sharp discontinuities but rather smoothly vary 
as a function of the collision energy The detailed predictions of Gyulassy and 



Rischke [44| about a slowly burning quark-qluon plasma have obtained a lot of 
experimental attention. Namely, it was expected that the HBT radius in the out 
direction will exceed the HBT radius in the side direction, predicting R ou t/ ' Rside > 1 
values reaching up to 4 or even 20 in case of an idea! b order phase —nfl 
In contrast, as is seen from Fig|2U R ou t ~ Rside ~ Riong in the RHIC energy region. 

The correlation function parameters may be related to physical ones characterizing 
the system size, lifetime, freeze-out duration, and expansion time. In particular, 

T 



based on a hydrodynamic expansion model it was shown that Ri ong = r / y ^ |i 
where Tf and T are the freeze-out time and temperature of the particles, respectively. 

One should emphasize that all these investigations suppose that centrality of heavy- 
ion collisions is under control and centrality scanning of the characteristics under 
discussion is an indispensable condition. 

Measurements of these quantities at the Nuclotron energies should be considered as a 
necessary continuation of global efforts to establish the energy dependence of properties 
of hadron production and search for signals of a phase change in nuclear collisions. 

In this connection the following theoretical and experimental studies at JINR are con- 
sidered as perspective 

1) research into the hadron properties in hot and/or dense baryonic matter. A spectral 
function change is expected, first of all for the cx-meson as a chiral partner of pions, 
which characterizes a degree of chiral symmetry violation. The rare specific channels 
of p-meson decays are also quite attractive. 

Solving these issues assumes a proper understanding of reaction mechanisms of high- 
energy colliding ions, knowledge of properties of strongly interacting QCD matter and its 
equation of state. In this respect, more general researches are in order: 

2) analyzing multiparticle hadron interactions, targeted to the development of a new 
statistical treatment as well as codes for space-time evolution of heavy nuclei colli- 
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sions at high energies. Particular attention should be paid to signals of a new phase 
formation during this evolution; 

3) studying the system size, lifetime, freeze-out duration, expansion time in the HBT 
analysis, scanning in atomic number and energy; 

4) analyzing the energy and centrality dependencies of the pion, hadron resonance and 
strange particle multiplicities, and the ratio of their yields, together with the trans- 
verse momentum, including K~, K*— and 0-meson spectra as well as manifestation 
of baryon repulsion effects on hadron abundances; 

5) studying dileptons (electron and muon pairs) production to see in-medium modifi- 
cation of hadron properties at high baryon densities; 

6) studying angular correlations in the transverse plane as well as radial, directed and 
elliptic flows; 

7) analyzing fluctuations of multiplicities, electric charge, and transverse momenta for 
secondary particles (their energy dependencies could give information on the phase 
transition range); 

8) analyzing nuclear fragments characteristics versus centrality, universality of nuclear 
fragmentation; 

9) energy and atomic number scanning for all characteristics of central heavy nuclei 
collisions (this might allow one to obtain information on the equation of state of 
strongly interacting QCD matter in the transition region), difference between central 
collisions of light nuclei and peripheral heavy ion collisions. 

The JINR Nuclotron has a possibility to accelerate heavy ions (up to A >200) to the 
maximal energy of 5 AGeV in about a year. This gives a chance to address experimentally 
many recent problems within the next several years before the FAIR GSI accelerator comes 
into operation. The proposed research program at the Nuclotron may be considered as a 
pilot study preparing for subsequent detailed investigations at SIS- 100/300 [l| and as an 
integral part of the world scientific cooperation to study the energy dependence of hadron 
production properties in nuclear collisions. 
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